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Abstract 

Following the scheme developed by Misner and Sharp, we discuss 
the dynamics of gravitational collapse. For this purpose, an interior 
cylindrically symmetric spacetime is matched to an exterior charged 
static cylindrically symmetric spacetime using the Darmois matching 
conditions. Dynamical equations are obtained with matter dissipat- 
ing in the form of shear viscosity. The effect of charge and dissipative 
quantities over the cylindrical collapse are studied. Finally, we show 
that homogeneity in energy density and conformal flatness of space- 
time are necessary and sufficient for each other. 

Keywords: Gravitational collapse; Dissipation; Junction conditions; 
Dynamical equations; Weyl tensor. 



1 Introduction 



A significant problem in gravitation theory and relativistic astrophysics is to 
understand the final fate of an endless gravitational collapse. A massive star 
undergoes a continual gravitational collapse at the end of its life cycle. This 
happens when a star has exhausted its nuclear fuel that provided a balance 
against the internal pull of gravity. 
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The importance of gravitational collapse in relativistic astrophysics was 
realized with the pioneer work of Oppenheimer and Snyder pQ. They used 
general relativity to study the dynamical collapse of a homogenous spherical 
dust cloud under its own gravity. Joshi and Singh |2] explored the spherically 
symmetric collapse of an inhomogeneous dust cloud. They found the end 
state of the gravitational collapse as a black hole or a naked singularity 
depending upon the initial density distribution and the radius of massive 
body. One may, however, consider dust as somewhat unrealistic form of 
matter, especially towards the end stages of a collapse, when pressures should 
be important. Keeping this fact in mind, the gravitational collapse of a 
perfect fluid and more general forms of matter have been studied. Misner 
and Sharp [3] discussed the gravitational collapse by taking the spherically 
symmetric ideal fluid in the interior and the Schwarzschild spacetime in the 
exterior of a star. They provided a full account of the dynamical equations 
governing the adiabatic relativistic collapse. 

Darmois [I] presented junction conditions that joined two solutions of the 
Einstein field equations across the surfaces of discontinuity. Sharif and Ah- 
mad [5J discussed junction conditions between static exterior and non-static 
interior spacetime in the presence of a positive cosmo logical constant. They 
also investigated the effect of a positive cosmological constant on spherically 
symmetric collapse with a perfect fluid. It was concluded that a positive 
cosmological constant slows down the rate of collapse. The same authors [B] 
also worked on cylindrical collapse of two perfect fluids using high speed ap- 
proximation scheme and examined the effects of pressure on the high speed 
collapse for two possible cases. Kurita and Nakao [7J discussed the collapse 
of null dust in the cylindrically symmetric spacetime and found a naked 
singularity at the symmetric axis. 

Gravitational collapse is a highly dissipative process [H]- [Hj] whose effects 
cannot be ignored in the study of collapse. Chan [11] studied a realistic 
model for a radiating star which undergoes dissipation in the form of radial 
heat flow and shear viscosity. He concluded that shear viscosity would in- 
crease anisotropy of pressure and also plays an important role in the study 
of gravitational collapse. The assumption of shear free motion of the fluid 
[12]- p]|] is used to obtain exact solutions of the field equations but it is an 
unrealistic approach. 

Herrera and Santos [9] studied dynamical description of gravitational col- 
lapse in view of Misner and Sharp's formulation. Matter under consideration 
was distributed with spherical symmetry and energy loss in the form of heat 
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flow and radiation. Herrera et al. [TJ] formulated the set of equations with 
regularity and matching conditions for the static cylindrically symmetric dis- 
tribution of matter. They showed that any conformally flat cylindrically 
symmetric static source cannot be matched to the Levi-Civita spacetime by 
using Darmois junction conditions. One of the authors (LH) [16J discussed 
the inertia of heat and its role in the dynamics of dissipative collapse. Herrera 
et al. [U] also formulated the dynamical equations to include dissipation in 
the form of heat flow, radiation, shear and bulk viscosity and then coupled 
with causal transport equations. Recently, Sharif and Rehmat [18] extended 
this work for the plane symmetric gravitational collapse. 

Some literature indicates keen interest for the inclusion of an electromag- 
netic field to discuss gravitational collapse. Bekenstein [T9] generalized the 
Oppenheimer-Volkoff equations of hydrostatic equilibrium [20] and Misner- 
Sharp formulation for the dynamics of spherical gravitational collapse to 
the charged case. Nath et al. [21] explored gravitational collapse in the 
presence of electromagnetic field by using the junction conditions between 
quasi-spherical Szekeres spacetime in the interior and the charged Vaidya 
spacetime in the exterior region. They concluded that formation of a naked 
singularity was enhanced by an electromagnetic field. Sharif and Abbas [22] 
investigated the effect of an electromagnetic field on the spherically symmet- 
ric collapse with the perfect fluid in the presence of positive cosmological 
constant. 

In a recent paper, Di Prisco et al. [23] derived dynamical equations for the 
spherically symmetric collapse by including an electromagnetic field. They 
concluded that Coulomb repulsion might prevent the gravitational collapse 
of the sphere. They also found the effect of charge on the relation between 
the Weyl tensor and the inhomogeneity of energy density. In this paper, we 
study the dynamics of a charged cylindrically symmetric spacetime to see 
the effect of charge on the rate of gravitational collapse. 

The format of the paper is the following. In the next section, we describe 
the gravitational source and some physical quantities. The Einstein-Maxwell 
field equation are given in section 3 and junction conditions are derived in 
section 4. We formulate the dynamical equations in section 5 and the relation 
between the Weyl tensor and the density homogeneity is given in section 6. 
The last section sums up the main results of the paper. 
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2 Interior Matter Distribution and Some Phys- 
ical Quantities 

We consider a cylindrical surface with its motion described by a timelike 
three surface E, which divides 4D spacetime into interior M~ and exterior 
M + manifolds. We assume co- moving coordinates inside hypersurface E. 
The interior cylindrically symmetric metric is given by 

ds 2 _ = -A 2 dt 2 + B 2 dr 2 + C 2 (# 2 + dz 2 ), 

— oo < t < +00, ^ r < +00, 0^0^ 2n, —00 < z < +00 (1) 

where {x~^} = {t, r, 4>, z} (/i = 0,1,2,3) and A, B and C are functions of 
t and r. Matter under consideration is anisotropic fluid which undergoes 
dissipation in the form of shear viscosity. The energy-momentum tensor for 
such a fluid is defined as 

T a p = (fi+ P±)V a Vp + P ±9a p + (P r - P ± )XaX? - 2Wafi, (2) 

where /i, P r , P±, r], V a and Xa are the energy density, the radial pressure, 
the tangential pressure, the coefficient of shear viscosity, the four-velocity 
of the fluid and the unit four-vector along the radial direction respectively. 
These quantities satisfy 

V a V a = -l, X a Xa = l, X a V a = 0. (3) 
The shear tensor a a b is defined by 

<7aP = V(a;[3) + Ch{aVp) ~ ^6(^ a/ g + V a Vp), (4) 

where the acceleration a a and the expansion are given by 

Oa = W, S = V£. (5) 

Let it be mentioned here that the bulk viscosity does not appear explicitly 
as it has been absorbed in the form of radial and tangential pressures of the 
collapsing fluid. The four-velocity and the unit four-vector are given by 

V a = A-^tift, x a = B- X S?. (6) 
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From Eqs.fll]) and ([H]), the non-zero components of the shear tensor are 



2 1 

0- U = -^ B2(J : a 22 = C 33 = --^zCV (7) 



The shear scalar a is defined by [24 



where dot refers to differentiation with respect to t. Thus we have 

a a0 a a0 = 2a 2 . (9) 
Using Eqs.flH} and (jSJ), it follows that 

A' n i (b c\ 

* = T 6= a{b +2 c)' (10) 

where prime represents derivative with respect to r. 

The C-energy for the cylindrical symmetric spacetime is defined by [25J 

E=]-(l-r 2 V a rV a r), (11) 

o 

where 

Here p is the circumference radius, I is the specific length, r is the areal 
radius, £ a stands for Killing vectors of cylindrically symmetric spacetime and 
E represents the gravitational energy per specific length of the cylinder. Thus 
the specific energy of the cylinder with the contribution of electromagnetic 
field in the interior region can be written as [26] 

F'- l + C ( d2 C ' 2 \+ S2 (12) 
E ~8 + 2 \7?-& (12) 



5 



3 The Field Equations 

The Maxwell equations are given by 

Fa/3 = 4>/3,a ~ 4>a,/3, (13) 
= /W", (14) 

where F a/3 is the Maxwell field tensor, (fi a is the four-potential and J a is the 
four-current. We can write the electromagnetic energy- momentum tensor in 
the form 

E <# = ^ (*2Ffh ~ \F^F lS g a ^ . (15) 

It is assumed that the charge is at rest with respect to the co-moving coordi- 
nates, thus the magnetic field is zero. Consequently, the four potential and 
the four current will become 

<p a = <M°, J a = P V a , (16) 

where $ = r) is an arbitrary function and p = p(t, r) is the charge 
density. The charge conservation yields 

rr 

s (r) = 2tt / pBC 2 dr, (17) 
Jo 

where s(r) is the total electric charge of the interior. For the interior space- 
time, using Eqs.flE]) and ([IB"]) , the Maxwell equations take the following form 



B' 


-2°- 


~B 


c 


B 


-2*' 


' B 


c 



>' = HopAB', (18) 

*'-(i + i- 2 cj*' = °- (19) 

Solving these equations simultaneously, it follows that 

The Einstein field equations for the interior metric can be written 

G aP = Bir(T a0 + E aP ). (21) 



6 



Using Eqs.Q, (j2J), (TT5]) and (|2T)j) . we can write the nonvanishing compo- 
nents as follows 

2 2 /l 2 

87r(T„„ + E„„) = 87r^.4 2 + *' 



4tt 2 C 4 

C \ B C J \B J \ C ' C V B C, 

(22) 

S,(T m + E m ) = = -2^--_ -— j, (23) 
8tt(Tii + En) = 8?r P r - -^r/a fi 2 Po 



BV/ 2^ + ^V-2^U^V + 2^' 



A/ l C \ C } AC I \C J AC 1 

(24) 



8n(T 22 + E 22 ) = 8ir[P ± + -^V°) C 2 + 



V3 "J 4tt 2 C 



2 



C\ 2 l B C A I B C\ BC 



A \B C A\B C BC 



c 








v c 


A 1 




i) 


C" 


- A '\ 






c 


A ' 




c 



C\* fA" C" A' ( B' C'\ B'C 



+ [b) {t + ^-a{b-c)- Sc)- (25) 

Equation (|23|) can be re-written in the following form 

1 C 

-(Q - V3a)' - V3a— = 0. (26) 

4 Junction Conditions 

In this section, we formulate the junction conditions for the interior and 
exterior manifolds. The interior manifold is given by Eq. (JT|) and the exterior 
manifold is the charged static cylindrically symmetric spacetime given by [27] 

ds\ = -NdT 2 + -^dB 2 + R 2 (d0 2 + dz 2 ), (27) 
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where 



R 



and x +M = {T, R, (p, z}. We can write the metric for the hypersurface £ in 
the following form 



(ds 2 



-dr 2 + A 2 (T)(d<p 2 + dz 2 



(28) 



where = (r, 0, z) (i = 0, 2, 3) represent the intrinsic coordinates of S. 
The Darmois junction conditions jl] can be stated as follows: 

• The continuity of the first fundamental form. This implies the conti- 
nuity of the metrics over the hypersurface 



(ds 2 )x = (ds 2 _h = (d4) E . 



(29) 



The continuity of the second fundamental form. This gives the conti- 
nuity of the extrinsic curvature over the hypersurface 



(30) 



Ktf is the extrinsic curvature given by 



, ^ d X id X ^ (ffjAi)I/ = 0jlj2>3)i (31) 



where are the components of outward unit normal to the hypersurface in 
the coordinates x ±fl - 

We can write the equations of the hypersurface as follows: 



/-(*,r) 
U(T,R) 



r — rs = 0, 

R — Re(T) = 0, 



(32) 
(33) 



where r s is a constant. Using Eqs. (l3"2"|) and (I3"3"i) in Eqs.Q and (12"7|) . we have 
the interior and exterior spacetimes on £ respectively 



(ds 2 _h 
(ds 2 + h 



A 2 (t, r^)dt 2 + C 2 {t, r E )(# 2 + <fe 2 ). 

/ / 1 i ■: . . \ 

N(Rv) - (N(Rx))- 1 



(34) 



dRy 



dT 



dT 2 + Ri(d<p 2 + dz 2 ). 



(35) 
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The continuity of the first fundamental form gives 



dt 
dr 

dT 
dr 



C(t,rv), 
1 

T 



dRj: 
dT 



(36) 
(37) 

(38) 



Now we consider the second fundamental form over E. For this purpose, we 
need the outward unit normals to S using Eqs. fl32l) and f|33l) . 



n, 



B(U, 1,0,0), 



(39) 

-f.1,0.0). (40) 



The surviving components of the extrinsic curvature K^j can be given as 



follows 



K, 



oo 



^o + o 



A' 
~AB 

d 2 TdR d 2 RdT N dN / dT\ 3 3 dN ( dR\ 2 dT 



K. 



22 



^2 + 2 



dr 2 dr dr 2 dr 




fcc\ 




^33 = 


\ B A 






dT' 




^3 + 3 = 


NR— 

dr 


s 



(41) 

(42) 
s 

(43) 
(44) 



The continuity of the second fundamental form, using Eq. (j30|) . yields 

3 O tat / m\ 2 

dr 2 dr dr 2 dr 2 dR \ dr J ' 2N dR \dr J dr 

A' 



d 2 TdR d 2 RdT NdNfdTy^ 3 dN (dR\ 2 dT 



dT' 
NR— 

dr 



AB 

CO 



(45) 
(46) 



Making use of Eq. ( 155|) , it follows 



d T 1 Lr fdR\ 2 

Tr- N\] N+ UJ' (47) 
Substituting Eq. fj4"T|) in Eq. (H6l) . we obtain 

Thus, if the interior and the exterior charges are equal over the hypersurface, 
i.e., Q = s, then we have 

E'- l - = M. (49) 

The difference between these two masses is equal to |, which is due to the 
least unsatisfactory definition of Thorne C-energy. Differentiating Eq. fj47|) 
with respect to r, then inserting this value in Eq. (I45j) and making use of 
Eq. fjloT) . we can write Eq.fT4"5l) as 

& BC A'C , . 

C-BC~AC =0 (50) 
This equation identically satisfies Eq. (l23j) . For the smooth matching of the 
interior and exterior metrics on hypersurface, Eqs.(l36 i) -( l38l) . ( f48i) and ( l50i) 
are the necessary and sufficient conditions. 

5 Dynamical Equations 

The conservation of energy- momentum, (T a P + E a ^).p = 0, implies that 

2C ( 2 \ 

ti + P± + - / =W 1=0 (51) 



and 



2C" / r \ /ihs' , , 
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In view of Misner and Sharp's formalism, we discuss the dynamics of a col- 
lapsing system. We introduce the proper time derivative as 



T Adt 



(53) 



The proper radial derivative Dr constructed from the circumference radius 
of a cylinder inside S is 

1 d 

Dr = ~ (54) 



R dr 



where 



R = C. (55) 
The fluid velocity in the case of collapse can be defined as 

U = D T (R) = D T (C) < 0, (56) 



which must be negative in the process of collapse. Using Eq. fjobj) . we can 
re- write Eq. ffl2|) as 



- 2 2 



If- + ^_ _ 1. I E' - - 



B ' C 2 C 

From Eqs. (j5^|) and fl55l) . Eq. fT26l) can be written as 



1/2 



(57) 



BE 



l ~D R {Q - V3a) - 



0. 



In non-dissipative shear free case, i.e., r] = a = 0, this equation takes the 
form 



0. 



This means that U ~ R which describes the homologous collapse. In view of 
Eqs. (|T2|) . (|23|) . (1241) and (|53|) . the rate of variation of the C-energy turns out 
to be 



D T E' = -AttR 2 ( P r - -j=pa - 



1 

32nR 2 



s 2 u ( $ i 



R 2 \8n 2 2 



(58) 
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The first term on the right-hand side of Eq. (158 j) in the case of collapse (U < 0) 
will increase the energy of the cylinder if 



(59) 



i.e., the effective radial pressure is greater than the particular value. This 
increase of C-energy is due to the work being done by the effective radial 
pressure. The second term in the round brackets describe energy leaving the 
system due to the Coulomb repulsive force. Similarly, using Eqs.f lT2|) . fl22|) . 
023]) and dHJ), we obtain 



D R E' = AnfiR 2 + - + —D R s + — 



8 R 



f4 



R 2 V 8tt 2 



(60) 



This equation indicates variation of the total energy between adjoining 
cylindrical surfaces inside the fluid distribution. The first term on right 
hand side gives the contribution of the energy density of the fluid element 
and the remaining terms are a constant and the electromagnetic contribution 
respectively. Integration of Eq.(l60l) leads to 



E' 



R 

4irfiR 2 dR + — + -— 



Mo 



8 2R 8vr 2 



R s 2 

w dR - 



(61) 



Using Eqs.( IT2|) . ( f24|) . ( )56|) and (!57|) . we can obtain the acceleration DtU of 
a collapsing matter inside Sas follows: 



D T U = - 



R 2 



E' 



I 



4irR(P r -— W ) 



EN 



/'o 



AB R 3 \8tt 2 



(62) 



Inserting the value of 4^ from Eq. fl62|) into Eq. fl52|) . it follows that 



fi + P r - -j=rio ] D T U 



1 

R 2 



E' 



I 



4tt R 



/' 
4 

7!' 



p ' ~ 7S" a 



rja J R — 



A + 1 

R 3 V 8vr 2 2 



-E 



Di 



P r - -^1°) + 2 (P r - P± - 2V3r]a) - 



l4s 



R lQir 3 R 4 



D R s 



(63) 
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This equation has the "Newtonian" form, i.e., 



Mass density x Acceleration = Force. 

Now we analyse the terms appearing in this equation as follows: The term 
in the round brackets on the left hand side represents inertial mass density. 
This gives the effect of the dissipative terms but there is no contribution of 
the electric charge. The remaining term on the left hand side is acceleration. 
There are two main terms on the right hand side. The first term represents 
the gravitational force. The factor within the round brackets is the same as 
on the left. It represents passive gravitational mass density by equivalence 
principle. The factor within the first square brackets shows how specific 
length, dissipation, and the electric charge affect the active gravitational 
mass term. Making use of Eq.flBTJ in Eq. (l6"3j) . it turns out 

R s 2 s 2 
W dR > R' 



which on differentiation yields 



| > D R s. (64) 



Thus the charge will increase the active gravitational mass only if this con- 
dition is satisfied. This increase of active gravitational mass causes the rapid 
collapse. 

The second term in the second square brackets constitute hydrodynamical 
forces. It consists of further three terms. The first contribution simply 
represents the gradient of the total effective radial pressure (including the 
influence of shear viscosity on P r ) which is always negative and is directed 
outward, this would prevent the gravitational collapse. The second term 
exhibits the effect of the local anisotropy of pressure with shear viscosity. If 
anisotropic pressure is positive then it increases the rate of collapse otherwise 
it decreases. In the last term, we have the Coulomb repulsion which may 
prevent the gravitational collapse of the cylinder. 

For hydrostatic equilibrium, i.e., U = 0, n — 0, Eq. (16"3"j) turns out to be 



2 

D R P r = - D R s - 2(P r - PA 



12 



R 2 V 8 7 R 3 V 8tt 2 2 



.(65) 
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We can also write Eq. (152"|) for the hydrostatic equilibrium U = 0, rj = as 
follows: 

* + > + + -5^ = 0. (66) 

This corresponds to the hydrostatic equilibrium for the spherically symmetric 
case and also gives the generalization of the Tolman-Oppenheimer-Volkoff 
equation for anisotropic charged fluid [28J. 

The static fluid leads to charged dust by taking P r = = Pj_, thus, it 
follows from Eqs. fl66|) and (1171) that 



A' ufapB 

»A-8^ = °- (67) 

Here B and C are function of r and hence B = C for a suitable transformation 
of r. Eliminating s from the field equations ( 1241) and ( 1231) . we get 

Air 2 

B = C, AB = l, S 2 = ^- T B' 2 . (68) 
Inserting these values in in Eq. (167j) . we obtain 

f = K P \ « = -£L. (69) 

We would like to mention here that such type of result was also found by 
Bonnor [29] for arbitrary spacetime. 



6 The Weyl Tensor 

Here we shall explore the relation between the Weyl tensor and density in- 
homogeneity. For this purpose, we define the Weyl scalar C 2 in terms of the 
Kretchman scalar 1Z, the Ricci tensor R a p and the curvature scalar R, i.e., 

C 2 = n-2R al3 R a ^ + ^R 2 . (70) 

Inserting the value of 1Z from Eq.( l88p in the appendix and making use of 
Eqs.p2>(123]) in Eq. (170]) . it follows that 
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where e is defined as 

e = -^R 3 . (72) 
485 

Applying the definitions of Dt and D R from Eqs. (I58|) and ( 1601) respectively 
to ( 17TT) . we obtain 

D T e = -47r[-R 3 D T ^-P r + P ± + 2V3ria S j 



and 



= 4n[~R 3 D R (ft - P r + P ± + 2V3 V a^ 



+ ( ft _P ± _ 2 ^)^] + £if-^ + ( s 



R \R 



(74) 



This shows that production of density inhomogeneity is directly linked to 
dissipative variables and charge. For the case of zero charge and dissipation, 
we have 

D R e + -^-R 3 D R f, = 0. (75) 

This implies that if D R fi = then C = (using the regular axis condition). 
Conversely, the conformally flat condition implies homogeneity in the en- 
ergy density. This result has already been verified for spherically symmetric 
gravitational collapse [23]. 



7 Summary and Conclusion 

To investigate how the system gradually changes with time, we have formu- 
lated a dynamical description of the cylindrically symmetric spacetime using 
Misner and Sharp's approach. Dissipative effects and anisotropic pressure 
have been taken into account. The junction conditions between cylindrically 
symmetric in the interior and charged static cylindrically symmetric space- 
time in the exterior provides the gravitational mass which causes gravity in 
the exterior region. 
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We have found the behavior of charge and pressure through the dynamical 
equations. It turns out that electric charge (unlike pressure) does not always 
produce a regeneration effect, i.e., the pressure trying to keep the star in 
equilibrium through the pressure gradients, at the same time contributes to 
the active gravitational mass. Thus it increases the gravitational attraction 
and hence it promotes stellar collapse at the same time. This is due to the 
inequality (63) which implies that if 

1 > DrS ' 

then it decreases the gravitational mass and also due to the Coulomb force 
that always opposes the gravitational force. We would like to mention here 
that our results indicate similarity with those found for the spherically sym- 
metric spacetime [23] . 

We have also established an expression indicating the relevance of the 
electric charge with the Weyl tensor and density inhomogeneity. Using the 
regular axis condition, it has been shown that homogeneity in energy den- 
sity and conformal flatness of spacetime are necessary and sufficient for each 
other. We would like to mention here that the Weyl tensor contains tidal 
forces that make the fluid more inhomogeneous in the process of evolution. 
It would be interesting to include also a heat flux and examine the corre- 
sponding transport equations. Also, one would be interested to extend these 
results for charged plane symmetric spacetime [30] . 

Appendix 

The interior metric has the following nonvanishing components of the Rie- 
mann tensor 

i?oioi = AA"-BB-^A'B' + jAB, (76) 

^0202 = -^(-CAB 2 + ACB 2 + AC'A 2 ), (77) 

#0212 = ^ {-C'AB + CA'B + BC'Aj , (78) 

#1212 = -^(-C"A 2 B + CBB 2 + B'C'A 2 ), (79) 

R2323 = ^(C 2 B 2 -C' 2 A 2 ), (80) 
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and 



R 



0202 



R 



0303 ) 



R 



0212 



r 



0313; 



R 



1212 



R 



1313- 



(81) 



Thus it has five independent components. The Kretchman scalar 1Z 

R"^ 5 R a/ 3-yS becomes 



TZ = A[ 



1 o 2 

(-Roioi) + 7777777(^0202)' 



(AS) 



(AG) 4 ' u u ; (AB) 2 G 4 
+ 7777777 (#1212) 2 + 777(^2323 



0212 j 



(82) 



The components of the Riemann tensor in terms of the Einstein tensor and 
the C-energy function can be written as 



-R0101 

-R0202 

-R0212 
-R1212 

-R2323 



(AB) 2 

(AC) 2 
C 2 



2A 2 



a 



00 



G11 



2B 2 G 3 



1 p , ^_22 2_ 

1 I 

8 2C 



E' 



I 



s 

2C 



-G 



01; 



2 

(5G) 2 



G, 



00 



1 



2A 2 G 3 
I 



E' 



8 2G 



(83) 
(84) 

(85) 
(86) 

(87) 



Inserting Eqs.(|83J-(|83) into Eq.(J82]), we get 
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I 



.2 \ 2 



16 



R C 6 \ 8 2Gj cA E ' 8 2G 



G, 



00 



A 2 



G\\ , G 



S 2 



22 



G 2 



G, 



01 



AS 



G 



00 



GopGn 
(AS) 2 + 



G11 

A 2 J ' ^B 2 
Gqo G11 \ G22 



+ 4 



G 



22 



G 2 



A 2 



S 2 



G 2 
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